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Abstract 

We present a comparative study of DNA nucleobases [guanine (G), adenine (A), thymine (T), and cyto- 
sine (C)] adsorbed on hexagonal boron nitride (A-BN) sheet and graphene, using local, semilocal, and van 
der Waals (vdW) energy-corrected density-functional theory (DFT) calculations. Intriguingly, despite the 
very different electronic properties of BN sheet and graphene, we find rather similar binding energies for the 
various nucleobase molecules when adsorbed on the two types of sheets. The calculated binding energies 
of the four nucleobases using the local, semilocal, and DFT+vdW schemes are in the range of 0.54 ~ 0.75 
eV, 0.06 ~ 0.15 eV, and 0.93 ~ 1.18 eV, respectively. In particular, the DFT+vdW scheme predicts not 
only a binding energy predominantly determined by vdW interactions between the base molecules and their 
substrates decreasing in the order of G>A>T>C, but also a very weak hybridization between the molecular 
levels of the nucleobases and the 7t-states of the BN sheet or graphene. This physisorption of G, A, T, and 
C on the BN sheet (graphene) induces a small interfacial dipole, giving rise to an energy shift in the work 
function by 0.11 (0.22), 0.09 (0.15), -0.05 (0.01), and 0.06 (0.13) eV, respectively. 

PACS numbers: 68.43.Bc, 82.39.Pj, 68.43. -h 
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I. INTRODUCTION 



The interaction of DNA nucleobases with inert surfaces has attracted much attention because 
of its importance for molecular recognition and self-organization processes.-"- Indeed, a great 
number of theoretical studies for the adsorption of the four DNA nucleobases [guanine (G), ade- 
nine (A), thymine (T), and cytosine (C)] on graphene have been performed to explore the binding 
mechanism and the relative binding strength of G, A, T, and C.— "— These essentially planar base- 
graphene model systems can simplify the structural complexity of a full three-dimensional DNA 
double-stranded or single- stranded polymer adsorbed on graphene, to provide a possibility for di- 
rect experimental characterization of the molecular interactions with graphene. According to pre- 
vious density-functional theory (DFT) calculations carried out within the local density approxima- 
tion (LDA),— the binding energy of the nucleobases on graphene varies in the following hierarchy: 
G>A?«T»iC. Calculations using the more accurate second-order M0ller-Plesset perturbation the- 
ory (MP2),— as well as other calculations utilizing DFT methods including van der Waals (vdW) 
interactions , 11112 obtained binding energy strengths of nucleobases with graphene in the ordering 
of G>A>T>C, consistent with the single solute adsorption isotherm study at the graphite- water 
interface. 16 The same trend was also confirmed using isothermal titration calorimetry.— In con- 
trast to the relatively large number of studies involving graphene, there have been relatively few 
studies concentrating on the interaction of DNA with the heterogeneous boron nitride (BN) sheet, 
which has been successfully fabricated through the micromechanical cleavage method 18 and the 
chemical- solution-derived method.— 

The hexagonal BN sheet exhibits the same honeycomb lattice structure as graphene. However, 
the electronic properties of the two sheets are drastically different from each other: graphene is a 
gapless semimetal with the nonpolar nature of the homonuclear C— C bond, while the BN sheet is 
an insulator with the polar nature underlying the charge transfer between its constituent B and N 
atoms.— ^ It is thus interesting to investigate and compare how such different electronic properties 
of the BN sheet and graphene affect the binding mechanism and the relative binding strength of the 
four nucleobases. Using LDA, Lin et al. showed that the nucleobase molecules are bound to the 
BN sheet via a polar electrostatic interaction, indicating that the interaction in the base-BN systems 
is somewhat different from the n-n interaction between the nucleobases and the nonpolar graphene 
sheet. 22 However, calculations based on LDA may not be reliable for this base-BN systems, since 
LDA cannot correctly describe the long -ranged vdW interactions. 
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In the present work, we investigate the adsorption of the four DNA nucleobases on the BN 
sheet and graphene, using vdW energy-corrected DFT (DFT+vdW scheme 23 ) calculations and, for 
comparison, LDA as well as generalized gradient approximation (GGA) calculations. Somewhat 
surprisingly, we find that the binding energy for a given base molecule is very similar on both the 
BN sheet and graphene. According to our analysis, it is revealed that despite the large differences 
in the individual atomic polarizabilities between the BN sheet and graphene, the vdW energy 
between adsorbed molecule and the BN sheet is nearly equal to the corresponding one in the 
base-graphene system, resulting in very similar contributions of vdW interactions to the binding 
energy. We find that the vdW interactions between base molecule and its substrate determine the 
sequence of the binding energy as G>A>T>C. Here, the contribution of vdW interactions to 
the binding energy amounts to ~1 eV, indicating a strong physisorption. Our electronic- structure 
analysis shows that this physisorption induces an interfacial dipole between the base molecule and 
substrate due to charge rearrangement, thereby causing a work-function shift relative to isolated 
BN sheet or graphene. 

II. COMPUTATIONAL METHODS 

The DFT+vdW,— GGA, and LDA calculations were performed using the FHI-aims 24 code 
for an accurate, all-electron description based on numeric atom-centered orbitals, with "tight" 
computational settings and accurate tier 2 basis sets. Calculations with larger tier 3 basis sets 
show that our binding energies are converged to within 0.01 eV. For the exchange-correlation 
energy, we employed the LDA functional of Ceperley- Alder— and the GGA functional of Perdew- 
Burke-Ernzerhof (PBE).— The DFT+vdW scheme was combined with the PBE functional, since 
PBE tends to underestimate the binding energy at organic/graphene interfaces. 27 In this PBE+vdW 
scheme, the total energy is composed of the PBE energy (£pbe) and the vdW energy (£Vdw) which 
is given by a sum of pairwise interatomic C(,R terms: 

£vdW = — ~ ^/damp(^AB,^^B^6AB^AB> (1) 
A,B 

where Rab is the distance between atoms A and B, C(,,ab is the corresponding C(, coefficient, R 1 ^ and 
i?g are the vdW radii, and /dampC^AB^A'^fi) is a damping function eliminating thei? A ^ singularity 
at small distances. Here, the C(, coefficients were computed in a first principles way from the PBE 
ground-state electron density using the recently developed DFT+vdW scheme.— 
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The present base-substrate systems were modeled in a periodic slab geometry in which base 
molecules adsorb on one side of the BN sheet or graphene. In order to make the overlap with 
electronic states from neighboring base molecules negligibly small, we employed a 5 x 5 unit cell 
of the BN sheet (graphene) with a vacuum spacing of 60 A between adjacent BN (graphene) 
sheets. The base molecules were terminated at the cut bond to the sugar ring with a methyl group 
to generate an electronic environment in the nucleobase more closely resembling the environment 
in an extended DNA chain rather than that of just individual isolated bases by themselves and also 
to introduce a certain degree of steric hindrance when interacting with the substrate. The k-space 
integration was done with a single T point in the Brillouin zone of the 5x5 unit cell. All atoms 
were allowed to relax along the calculated forces until every residual force component amounted 
to less than 0.02 eV/A. 

III. RESULTS AND DISCUSSION 

We begin by describing the equilibrium structures of adsorbed G, A, T, and C on graphene ob- 
tained using the LDA, PBE, and PBE+vdW schemes. On the basis of previous theoretical calcula- 
tions, -i&Iii 2 ^ we optimize the AB-stacking-like arrangements of all four nucleobases on graphene. 
Here, we use the optimized lattice constant of graphene as 2.445, 2.467, and 2.465 A (see Table 
I) for LDA, PBE, and PBE+vdW, respectively. The optimized AB-stacking-like structures for ad- 
sorbed G, A, T, and C on graphene are very similar to the corresponding ones on the BN sheet (see 
Fig. 1). Table II lists the vertical distance between the base molecule and graphene sheet, obtained 
using LDA, PBE, and PBE+vdW. We find that LDA, PBE, and PBE+vdW give the values of the 
vertical distance ranging 3.08~3.17 A, 3.95~4.02 A, and 3.26~3.29 A, respectively. Thus, the 
PBE+vdW values of the vertical distance are seen to fall between the LDA and PBE ones. 

The calculated binding energies for the four base-graphene systems using LDA, PBE, and 
PBE+vdW are listed in Table III. We find that LDA (PBE; PBE+vdW) gives E h = 0.72 (0.14; 
1.18), 0.55 (0.06; 1.00), 0.54 (0.08; 0.95), and 0.56 (0.13; 0.93) eV for adsorbed G, A, T, and 
C on graphene, respectively. Thus, the binding sequence of the four nucleobases is G>Af=aT^C 
for LDA and G>A>T>C for PBE+vdW. This trend of the binding sequence for the four base- 
graphene systems agree well with previous theoretical results (see Table III),— - — though the mag- 
nitudes of £b vary depending on the employed implementations of various computational meth- 
0( j s 28,29 j t j s no t a ble that PBE gives lower binding energies of less than ~0.2 eV (Table III) and 
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(a) (b) 




FIG. 1 : (Color online) Top and side views of the optimized structures of adsorbed (a) G, (b) A, (c) T, and 
(d) C on the BN sheet, respectively. The circles represent B, C, N, O, and H atoms with decreasing size. 
For distinction, N atoms on the BN sheet (nucleobases) are drawn in light (dark) color. The vertical plane 
along the dashed line is taken for the contour plot in Fig. 4. 

relatively larger vertical distances around ~4.0 A (Table II), indicating that it does not correctly 
describe the n — n stacking interactions between the base molecule and graphene. We also note 
that LDA predicts relatively lower binding energies yet smaller vertical distances compared with 
those obtained using PBE+vdW. 

To assess the contribution of vdW interactions to the binding energy, we decompose the 
PBE+vdW binding energy into two parts computed from £pbe and £VdW : 

Eb = £b,PBE + £b,vdW, (2) 

where £ b ,PBE=-[£pBE(base/sub) -E PB E(base) -E PB E(sub)] and £ bjVdW =-[£ vdW (base/sub) - 
Evdw(base) — E vc rw(sub)], in which "sub" stands for "substrate". Figure 2 shows £b,vdw for the 
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TABLE I: Optimized lattice constant (in A) of the BN and graphene sheets obtained using LDA, PBE, and 
PBE+vdW, in comparison with previous theoretical results. 



LDA 



PBE 



PBE+vdW 



BN 



graphene 



This 
Ref. 32 
Ref. 33 
This 
Ref. 30 
Ref. 31 



2.489 
2.49 

2.445 
2.445 



2.513 

2.511 
2.467 

2.47 



2.510 



2.465 



TABLE II: Calculated vertical distance (in A) between the base molecule and BN or graphene sheet, 



obtained using LDA, PBE, and PBE+vdW. 

G A T C 

BN LDA 3.03 3.04 3.08 3.12 

PBE 3.80 3.94 3.96 4.04 

PBE+vdW 3.21 3.25 3.27 3.24 

graphene LDA 3.08 3.17 3.10 3.12 

PBE 3.95 4.00 4.02 3.97 

PBE+vdW 3.26 3.29 3.29 3.27 



four base-graphene systems. We find that adsorbed G, A, T, and C on graphene have fi^vdw 
(£b,PBE) = 1.18, 1,10, 1,02, and 0.94 (0.00, -0.10, -0.07, and -0.01) eV, yielding E b = 1.18, 
1.00, 0.95, and 0.93 eV, respectively. Thus, for all four nucleobases, the contribution of vdW 
interactions to the binding energy amounts to ~1 eV, indicating strong physisorption. We note 
that the PBE+vdW binding energies for adsorbed G, A, T, and C on graphene are larger than 
the PBE ones by 1.04, 0.94, 0.87, and 0.80 eV, respectively (see Table III). These values of the 
binding energy difference between the PBE+vdW and PBE calculations are somewhat different 
from the corresponding ones of £b,vdw due to the use of two different adsorption structures in 
the PBE+vdW and PBE calculations (see Table II). It is likely that the relatively shorter vertical 
distance in PBE+vdW compared with that in PBE gives rise to the Pauli repulsion between n 
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electrons of the nucleobase molecule and graphene, yielding a negative value for £b,PBE- 




Q) 0.8- 
- 0.6- 
% 0.4- 

uf 21 



FIG. 2: (Color online) Calculated £b,vdw for adsorbed G, A, T, and C on the graphene (top panel) and BN 
sheets (bottom panel). M-C represents the value for the base-graphene systems. For the base-BN systems, 
the two components are decomposed: M-B (M-N) represents the component arising from the atomic pairs 
between the constituent atoms in the base molecule and the B (N) atom. 



TABLE III: Calculated binding energies (in eV) of G, A, T, and C adsorbed on the BN and graphene sheets, 
in comparison with previous theoretical results. 









G 


A 


T 


C 


BN 


LDA 


This 


0.75 


0.56 


0.57 


0.59 






Ref. 22 


0.69 


0.58 


0.56 


0.54 




PBE 


This 


0.15 


0.07 


0.08 


0.13 




PBE+vdW 


This 


1.18 


1.01 


0.94 


0.93 


graphene 


LDA 


This 


0.72 


0.55 


0.54 


0.56 






Ref. 10 


0.61 


0.49 


0.49 


0.49 




PBE 


This 


0.14 


0.06 


0.08 


0.13 




PBE+vdW 


This 


1.18 


1.00 


0.95 


0.93 






Ref. 12 


0.99 


0.85 


0.76 


0.76 




vdW-DF 


Ref. 12 


0.74 


0.63 


0.60 


0.58 



In Fig. 2, it is clearly seen that the binding energy sequence G>A>T>C is determined by 
vdW interactions between the base molecule and graphene. Since the vdW energy in the present 
PBE+vdW scheme is given by a sum of pairwise interatomic C(,R~ 6 terms, we can say that the 
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magnitude of the effective C^u coefficient between the base molecule and the C atom of graphene 
is in the same order of G>A>T>C as the binding energy, consistent with a previous theoretical 
study 10 that the strength of the binding energy is governed by the polarizabilities of the base 
molecules which are in the order of G>A>T>C. 

Next, we study the adsorption of four nucleobases on the BN sheet using the LDA, PBE, and 
PBE+vdW schemes. The optimized PBE+vdW structures for adsorbed G, A, T, and C on the BN 
sheet are respectively shown in the panels (a), (b), (c), and (d) of Fig. 1, showing the same AB- 
stacking-like arrangement as the case of the base-graphene systems. In Table II, we find that LDA, 
PBE, and PBE+vdW give the values of the vertical distance ranging 3.03~3.12 A, 3.80~4.04 A, 
and 3.21~3.27 A, respectively. These values in the base-BN systems are slightly smaller than 
the corresponding ones in the base-graphene systems (see Table II). The calculated LDA (PBE; 
PBE+vdW) binding energies of G, A, T, and C on the BN sheet are 0.75 (0.15; 1.18), 0.56 (0.07; 
1.01), 0.57 (0.08; 0.94), and 0.59 (0.13; 0.93) eV, respectively. It is seen that LDA predicts a 
slightly larger binding energy on the BN sheet, compared with that on graphene (see Table III). 
This may be attributed to a more electrostatic attraction between the nucleobase and BN sheet, due 
to the polar nature of the B— N bond. On the other hand, either PBE or PBE+vdW gives almost 
the same binding energy on the graphene and BN sheets, indicating that each base molecule binds 
to the two sheets with a nearly equal binding strength. Indeed, PBE+vdW predicts that adsorbed 
G, A, T, and C on the BN sheet have £b,vdw = 1-16, 1.09, 0.99, and 0.93 eV, respectively. These 
values are very close to the corresponding ones (1.18, 1,10, 1,02, and 0.94 eV) on graphene. To 
understand this similarity of £t>,vdw on the two substrates, we decompose £b,vdw obtained on the 
BN sheet into two components: i.e., one (the other) component arising from the atomic pairs 
between the constituent atoms in base molecule and the B (N) atom. 34 Such decompositions for 
the four base-BN systems are displayed in Fig. 2. We find that the ratio of magnitudes of the 
two components in all base-BN systems is about 8:5, indicating that the effective C(, : ij coefficient 
between the base molecule (M) and the B atom is greater than that between M and the N atom 
by a factor of ~1.6. Since the magnitudes of £b,vdw in the base-graphene and base-BN systems 
are very close to each other, the value of (C6.m-b + Q,m-n)/2 can be nearly equal to the value of 
Q,M-c- The interesting conclusion from these results is that, despite the different bonding natures 
(i.e., non-polar and polar) of the graphene and BN sheets, the vdW interactions between each 
nucleobase molecule and the two substrates are close to each other, resulting in similar binding 
energies in the base-graphene and base-BN systems (see Table III). It is noteworthy that the in- 
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plane polarizability of graphene and BN sheet would be very different, in the sense that the more 
localized electronic density of the polar network of B- and N-atom in the latter is more difficult 
to deform (from an energetic point of view) than in the more delocalized electronic density of 
the non-polar network of C-atoms in graphene. Despite this in-plane difference in polarizability, 
the main contribution to the vdW interaction with the tc— TC stacked nucleobases seems to arise 
however mainly from the eponymous 71-orbitals belonging to the substrate (which are much closer 
in distance to the corresponding ft-orbitals of the nucleobases). It is thus likely that the 71 — K 
stacking interactions on the BN sheet will be very close to those on graphene, resulting in the 
overall quantitatively similar binding energy for a given nucleobase on both the BN sheet and 
graphene. 

(a) 

4 
3 
2 

I 1 
D) 

is 

c 

LU -1 

-2 
-3 
-4 

K r M K K r M KK rMKK r M K 

(b) 
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3 
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D) 

is 

c 
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-2 
-3 
-4 

K TMKK r M KK r M KK TMK 

FIG. 3: (Color online) Calculated band structures of adsorbed G, A, T, and C on the (a) graphene and (b) 
BN sheets. The band dispersions are plotted along the symmetry lines of the Brillouin zone of the 5x5 unit 
cell [see the inset in (a)]. The energy zero represents the Fermi level. For the base-BN systems, the Fermi 
level lies in the middle of the band gap of the BN sheet. For distinction, the highest occupied molecular 
orbital (HOMO) and lowest unoccupied state (LUMO) are drawn with thick lines. 

Figure 3(a) and 3(b) show the calculated PBE+vdW band structures for adsorbed G, A, T, and 
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C on the graphene and BN sheets, respectively. It is seen that the molecular orbitals of the four 
nucleobases hardly hybridize with the 71 states of the graphene or BN sheet. We note that the 
HOMO of adsorbed G, A, T, and C on the graphene (BN) sheet locates at 0.93 (1.80), 1.23 (2.04), 
1.53 (2.31), and 1.33 (2.17) eV below the Fermi level, respectively. Here, the HOMO positions 
follow the order of the ionization energies for the isolated base molecules G, A, T, and C which are 
calculated to be 5.39, 5.48. 5.79, and 5.60 eV, respectively. Our analysis of the Mulliken charges 
in the PBE+vdW calculations shows a very small charge transfer of less than 0.03 e between any 
of the four nucleobases and the graphene or BN sheet. To see the rearrangement of charge at the 
base-substrate interface, we calculate the charge density difference defined as 

Ap = PM/Sub-(PM + Psub), (3) 

where PM/sut>> Pm> and p su b denote the charge densities of the base-substrate system and its sepa- 
rated systems, i.e., isolated layer of base molecules and clean substrate, respectively. Figure 4(a) 
and 4(b) show Ap for adsorbed G on the graphene and BN substrates, respectively. It is seen that 
the two base-substrate systems involve charge rearrangement at the interface. This charge rear- 
rangement gives rise to an interfacial dipole, thereby causing the work-function shift.— - — The 
calculated work-function shifts of all base-substrate systems are listed in Table IV. We find that, 
upon adsorption of G, A, T, and C on the graphene (BN) sheet, the work function increases relative 
to the value of 4.24 (3.48) eV at the isolated graphene (BN) sheet by 0.22 (0.1 1), 0.15 (0.09), 0.01 
(—0.05), and 0.13 (0.06) eV per molecule within a 5x5 unit cell, respectively. We note that the 
work function of isolated graphene is calculated to be 4.24 eV, which is somewhat underestimated 
compared to the experimental 38 value (~4.6 eV) for graphite but consistent with other PBE calcu- 
lations.— The work-function shift AW can be correlated with the induced interfacial dipole Ap by 
a simple electrostatics relation AW = eAp/(£oA)j 40 where A is the area of 5 x5 unit cell. Using this 
relation, we estimate Ap as 0.065 (0.032), 0.045 (0.028), 0.002 (-0.015), and 0.040 (0.019) ek 
for adsorbed G, A, T, and C on the graphene (BN) sheet, respectively. Note that Ap corresponds 
to the normal component of induced dipole moment directed from base molecule to substrate. On 
the basis of our PBE+vdW results, we can say that, although the binding mechanism between the 
nucleobases and the graphene or BN sheet is driven by the vdW interactions, the interfacial dipole 
is induced upon adsorption to yield the work-function shift in the order of G>A>C>T. Here, we 
note that the values of Ap for adsorbed T on the graphene and BN sheets are close to zero, possibly 
because of the cancelation of inhomogeneous interfacial dipole moments. 
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(a) 




FIG. 4: (Color online) Calculated charge density difference Ap for adsorbed G on the (a) graphene and (b) 
BN substrates. The contour plot in (a) is drawn in the vertical plane along the dashed line in Fig. 1(a). 
The same vertical plane is also used in (b). The first solid (dashed) line is at 0.0004 (—0.0004) el A 3 with 
spacings of 0.0004 e/A 3 . 



TABLE IV: Calculated work-function shift (in eV) upon physisorption of nucleobases G, A, T, and C on 
the BN and graphene sheets. The reference work functions of pristine BN and graphene sheets are given in 
the last column (in eV). 





G 


A 


T 


C 


pristine 


BN 


+0.11 


+0.09 


-0.05 


+0.06 


3.48 


graphene 


+0.22 


+0.15 


+0.01 


+0.13 


4.24 



IV. SUMMARY 



We have investigated the adsorption of the four DNA nucleobases on the BN sheet and on 
graphene, using the LDA, PBE, and PBE+vdW schemes. The calculated binding energies of the 
four nucleobases on the two different substrates were predicted in the order of PBE+vdW > LDA 
> PBE. We found that the vdW interactions between each base molecule and the two sheets are 
very close to each other, giving rise to similar binding energies in the base-BN and base-graphene 
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systems. Here, the magnitudes of the vdW interactions range from 0.9 to 1.2 eV, indicating a 
strong physisorption. We also found that the variation of vdW interactions depending on the base 
molecules determines the sequence of the binding energy as G>A>T>C, following the hierarchy 
of polarizabilities of the four DNA nucleobases. Our analysis shows that this physisorption induces 
an interfacial dipole between the base molecule and the substrate, leading to a small change in the 
work function relative to isolated graphene and BN sheets. 
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